Induced p-wave superfluidity in strongly interacting imbalanced Fermi gases 
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The induced interaction among the majority spin species, due to the presence of the minority 
species, is computed for the case of a population-imbalanced resonantly-interacting Fermi gas. It 
is shown that this interaction leads to an instability, at low temperatures, of the recently observed 
polaron Fermi liquid phase of strongly imbalanced Fermi gases to a p-wave superfluid state. We find 
that the associated transition temperature, while quite small in the weakly interacting BCS regime, 
is experimentally accessible in the strongly interacting unitary regime. 

PACS numbers: 05.30.Fk, 03.75.Ss, 67.85.-d, 32.30.Bv 



The extraordinary controllability of cold atomic gases 
has yielded a wide range of interesting phases of matter, 
including a bosonic Mott insulator and a paired super- 
fluid state of two species of atomic fcrmions [1, 2]. In the 
latter setting, experiments have demonstrated control of 
both the interspecies interactions and the relative density 
of the two spin states [3, 4], with the latter experimen- 
tal knob being deleterious to pairing and superfluidity, 
which favors an equal density of the two species. 

Thus, experiments on such imbalanced Fermi gases can 
probe the stability of superfluidity in a correlated system 
and therefore may shed light on the tendency towards 
pairing in related systems, such as electronic supercon- 
ductors. The phase diagram of imbalanced Fermi gases is 
quite rich [5, 6], possessing regions of imbalanced super- 
fluidity, phase separation, normal Fermi liquid, and (pos- 
sibly, though not yet observed) a region of exotic Fulde- 
Ferrell-Larkin-Ovchinnikov (FFLO) superfluidity [7]. 

Our present focus is on the strongly imbalanced region, 



when the polarization P 



(with n a the density of 



fermion species a) is close to unity, with a small density 
of spins- J, immersed in a spin-f Fermi sea. Experiments 
in this regime [8, 9] have found results consistent with 
the formation of spin polarons [10], in which a cloud of 
spins-! gather around each spin- J,, leading to a polaron 
Fermi liquid state [11-13]. 

The polaron theory of the strongly imbalanced Fermi 
gas predicts that both the spins- f and spins-! are Fermi 
liquids, exhibiting sharp Fermi surfaces at low temper- 
ature T. However, general arguments due to Kohn and 
Luttingcr (KL) [14, 15] predict that such Fermi surfaces 
must be unstable to other ordered states as T — > 0. 
A natural question then emerges: What is this ordered 
state for imbalanced Fermi gases? Since the imposed im- 
balance (and concomitant Fermi-surface mismatch) pre- 
cludes s-wave singlet interspecies pairing, one instead ex- 
pects triplet (likely p-wave) intraspecies pairing [14-16] 
of the spin-! and spin- 1 fermions. 

In this Rapid Communication we develop a theoretical 
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FIG. 1: (Color online) The p-wave channel of the effective 
interaction, determined by the diagrammatic inset, between 
majority spins, ?/ =1 (fcp^-, k-p-f), multiplied by the Fermi-energy 
density of states iV^(e F ^), as a function of the density im- 
balance P and s-wave scattering length a B , is shown. Here, 
k F — (k F ^ + fc F jJ/2. At large P, above the dashed line, it 
is attractive leading to a p-wave superfluid at temperatures 
below T c . The solid white line labels the location of a line 
of FFLO quantum critical points. This coincides with the 
location where v =1 (fc F t, &Ft) is large. 



description of the induced interactions among the ma- 
jority spin-! fermions in the presence of a small density 
of spins-!, over a broad range of interaction and polar- 
ization values. We find an attractive effective interac- 
tion at the spin-! Fermi surface, shown in Fig. 1, leading 
to an instability towards p-wave superfluidity below a 
transition temperature T c , computed below [17]. In the 
extreme weak coupling Bardcen-Coopcr-Schricffcr (BCS) 
limit, where the s-wave scattering length a s — > 0~, the 
p-wave superfluid transition temperature has been com- 
puted [16]; unfortunately, it is exceptionally small (in 
agreement with KL [14]), with T c oc exp[— c/(kF^-a 3 ) 2 }, 
where fcp-|- is the spin-! Fermi wave vector and c is a 
constant of order unity. The pairing mechanism is quite 
simple in this perturbativc limit: density fluctuations of 
one species leads to an attraction between particles of op- 



2 



posite spin. For spins-t, this induced interaction is pro- 
portional to the spin- 4. density-density correlation func- 
tion (i.e., the Lindhard function). Although a p-wave 
supcrfluid is predicted in the BCS limit, most experi- 
ments occur in the unitary region where the interspecies 
interactions are strong, \a s \ — > 00, invalidating simple 
perturbative results. Our analysis of induced interac- 
tions in the unitary regime involves extending the ladder 
approximation (known to describe the polaron Fermi liq- 
uid regime discussed above) to include subleading classes 
of Feynman diagrams. As seen in Fig. 1, the predicted 
effective interaction in the p-wave channel can be quite 
large near the unitary regime We find a maximum transi- 
tion temperature of fc B T c ~ 0.03e F ^., which is low but not 
unreasonable given recently reported temperature scales 
(e.g., Rcf. 18). 

Our starting point is the standard one-channel model 
for two species (er =1", 4-) of fermion (c^ ) interacting via 
an ,s-wave Feshbach resonance [19]. The Hamiltonian is 



A 



/ j ^kcr^kcrHcCT 
k : er 



V l -k+qt L k'-qi L k'l'-'kt! 
k,k',q 



(1) 



where £ kcr = £k — (i a , with dispersion — k 2 /2m 
(H = 1) and chemical potential fi a ; V is the system 
volume (henceforth set to unity), and A is the coupling 
strength of a short-ranged pseudo-potential, related to 



the experimentally controllable (via a magnetic field- 
tuned Feshbach resonance) scattering length a s via 



47ra„ 



= A" 



A 1 



(2) 



where A is an ultraviolet cutoff; below we shall take the 
limit A — > 00 while A — > _ , such that a s remains fixed. 

We are interested in the phases of Eq. (1) in the 
strongly imbalanced limit P — > 1 and proceed (in 
the spirit of mean-field theory) by assuming the pres- 
ence of a self-consistently determined pairing ampli- 
tude A-|-(k, (j) among the spins-f, but not the spins- 
4 [17]. Under this assumption, the 2x2 spin-f 
Nambu Green's function, in Fourier-Matsubara space, is 
'Gf(K) F t (K) 
F t (K) -G t (-K) y 

K = (zw„,k). Gf(K) satisfies the Dyson equation [20] 
Q^(K) = G^l(K) - t^(K), with bare Green's func- 

tion y tfi (K) = 1 



where the four-vector 







and self-energy 



/ E+(K) -A-t-(K) \ 
Ef(K)= I a* (if) -£ t (-A)J' Using the equation of 

motion [21], the self-energy can be expressed in terms of 
the two-particle vertex F, by 



fit (Jf) = Xa z Gi(Ki) [a - ^ g^{K 2 )G x {K + Kx- K 2 )f(K 2 , K + K x - K 2 , K\,K) 

K x K 2 



(3) 



r 



where <jq is the 2x2 identity matrix, <r z is a Pauli matrix, 
and G±(K) is the scalar spin-4 Green's function, which 
satisfies a similar expression, but with a diagonal self- 
energy (since we assumed the spins-4 are unpaired) . The 
summation is short for Y] K = k B T^2 iuj Y]^,. 



Although Eq. (3) is in principle exact, to make progress 
we must make a physically motivated approximation for 
r, corrcpsonding to certain classes of Feynman diagrams. 
Previous work has analyzed the phases of imbalanced 
Fermi gases within the i-matrix or ladder approxima- 
tion [11, 22, 23]. The set of diagrams associated with 
the ladder approximation also emerges in the large- N ap- 
proximation [12], in which one generalizes the model to 
consist of 2N species of fermion. Within the present for- 
malism in which Q^(K) has Nambu structure, these con- 
tributions arise from including the ladder plus the maxi- 
mally crossed self-energy diagrams, sketched in Figs. 2(a) 
and 2(b), respectively. We first analyze Eq. (3) including 
only these diagrams. Exchanging A for a s using Eq. (2), 
and keeping only contributions that arc finite in the limit 



A — > 00, we find a diagonal self-energy; 







-t(Q - K] 



(4) 



indicating the absence, at this level, of pairing for the 
spins-f. Here, 



t(K)- 1 = 



^G^-W)^! (5) 



is the usual scalar i-matrix. 

Thus, the contributions from Figs. 2(a) and 2(b) yield 
an unpaired solution for the self-energy, i.e., a Fermi 
liquid, as found by previous i-matrix or large-A theo- 
ries [11, 12, 22, 23]. Given that the KL arguments im- 
ply the eventual instability of this state, we now turn 
to subleading contributions to the self-energy, shown in 
Fig. 2(c), that possess ladder and crossed-laddcr subdi- 
agrams. Again replacing A for a s using Eq. (2), we find 
that these diagrams yield an off-diagonal contribution to 
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FIG. 2: The ladder (a) and maximally crossed (b) diagrams 
contributing to the self-energy, which define the ladder t and 
crossed-ladder t c diagrams. The bottom row, (c), shows a 
subleading class of diagrams with both ladder and crossed- 
ladder diagrams that lead to a pairing instability. The spin-f 
lines are 2x2 Nambu Green's functions Q^(K), spin- 4. lines 
are normal scalar Green's functions G±(K), and the wavy 
lines are Act.. 

the self-energy, i.e., a pairing amplitude A^(K) given by 

A t (K)=Y / V(K,K')F t (K') 1 (6) 
Q 

V(K, K') =]T t[P-K)tip+K')G^G^K'+P -K){7) 
p 

Here, V(K, K') is the effective induced interaction among 
spins-t; the corresponding Fcynman diagram is shown 
in Fig. 1. In principle, the integral equation (6) must 
be solved sclf-consistently along with the diagonal self- 
energy, Eq. (4). However, we shall we use some physi- 
cally motivated approximations to simplify our analysis, 
focusing on the onset of pairing of the spins-t at a tem- 
perature T c (above which the system is a Fermi liquid). 
We assume the presence of a static momentum-dependent 
pairing order parameter, A^(K) = A^.(k), and neglect 
the frequency dependence of V(K,K') [15]. For the di- 
agonal components of the self-energy, we simply assume 
that the chemical potential is renormalized to the Fermi 
energy via ji a —t \i a — S^k^O) = e-p a (consistent with 
the Luttinger theorem [24, 25]). Within these approxi- 
mations and after analytic continuation the effective in- 
teraction takes the form 

V(k, k') « 2Rc[ ^'(k + q, £qi)i a (q - k', U) (8) 
q 

xGi(k-k' + q,£ qi )n F (£ qi )], 

where r/a refers to the retarded or advanced quantities 
and n F (uj) is the Fermi distribution function. Equation 
(6) then simplifies to 



A t (k) = -^y(k,k')^^tanh 



2T 



with i?k = y^kt + l^tWI 2 ' the solution of which re- 
quires an understanding of the momentum structure of 
the effective interaction V(k, k') in the vicinity of the 
spin-f Fermi surface. The transition temperature T c for 
a given angular momentum is found by solving the lin- 
earized, in A.j.(k), version of Eq. (9) and projecting onto 
the relevant channel [19]. Assuming p-wave pairing, one 
needs the £ = 1 projection of the induced interaction, 
v e=1 (k 7 k') = Jq d9 sin cos ^^(k, k'), where is the an- 
gle between k and k'. Furthermore, we find via a direct 
numerical integration of Eq. (8), that v 1 (k, k') is only ap- 
preciable for k and k' within a range k-p\, of each other; 
this defines an effective bandwidth, of the order of cfj,, 
over which the induced interaction is nonzero. 

The strong induced attraction among the spins-f, 
shown in Figs. 1 and 3a, suggest a robust p-wave super- 
fluid at T — > 0; to estimate the associated transition tem- 
perature T c we must make further approximations. The 
remaining momentum integrations in Eq. (9) arc sharply 
peaked at the Fermi surface, yielding the result 



kv>T r 



2cT 



-e P | exp 



iV t (e F t)^ _1 ( fc Ft ; fc Ft) 



(10) 



with 7 the Euler gamma constant. We have also 
found [26] the same result for the transition tempera- 
ture (and the same effective interaction, Eq. (7)) via a 
somewhat different approach by considering the Thou- 
less criterion [27] for T c , determined by the point at which 
the spin-f pair-pair fluctuations in the normal state be- 
come unbounded. Within such an approach, Eq. (7) is 
the irreducible vertex in the particle-particle channel of 
the Bcthc-Salpeter equation [28]. 

In general, (8) has to be determined numerically, but 
analytic results can be found for certain limiting cases. 
In the asymptotic BCS limit a s — > 0~ , the i-matrix 
£ r / a (k, lj) — > 47ra s /m, and Eq. (8) reduces to the result 
of Ref. [16] (where, as we have noted, T c is extremely 
small). Analytic results can also be obtained in the ex- 
tremely imbalanccd limit, i.e., kp^/k-p^ 3> 1. In this limit 
the i-matrices appearing in Eq. (8) also become indepen- 
dent of k and k'; t^(k n ± q,£ qi ) -+ - ■ 

Evaluating the remaining integral gives 



2e 7 

k B T c « eF+exp 

7T 



2 \n(z) \ 2k F ^a 



2 



(11) 



(9) 



with z = k-pf/kpi. This formula correctly captures the 
vanishing of T c for P — > 1, but doesn't adequately cap- 
ture the peaks shown in Fig. 3, which where found by a 
direct numerical analysis of Eq. (8). These results indi- 
cate the presence of pairing at an experimentally acces- 
sible temperature in unitary imbalanced gases. 

The peak in the induced attraction, and in T c , at large 
P, can be understood by noting that a crucial contri- 
bution comes from particle-hole excitations at the spin- J, 
Fermi surface. However, particlc-holc excitations with 
a transferred momentum larger than 2fcpj. are energeti- 
cally suppressed, so that the associated density response 
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FIG. 3: The top panel shows the 1 = 1 channel of the effective 
interaction for spin-t fermions, times the density of states, as 
a function of polarization, at unitary and into the BEC side. 
The bottom panel shows the corresponding p-wave transition 
temperature, according to Eq. (10), with eFt/eFi = [(1 + 
P)/(l — P)] 2//3 , via a numerical integration of Eq. (8). 

exhibits a strong momentum dependence at k ~ 2fcp.|., 
and a correspondingly large p-wave projected interaction 

k 3 — k 3 

at fcpf ~ or?= , j) T ,i l ~ 0.78. Near unitarity, 

K Ft" 1 " F4. 

the magnitude of this peak is enhanced by the presence of 
strong fluctuations toward the FFLO phase, signaled by 
a divergence of the t- matrix at a non-zero wave vector Q, 



i.e., i r (Q,0) — > oo. To illustrate this, in Fig. 1 the white 
line shows the P at which a quantum phase transition to 
the FFLO state occurs; thus, below this line, the p-wave 
paired phase may undergo a second phase transition to 
the FFLO. 

The confirmation of our scenario will require detecting 
the onset of p-wave pairing at T c and the properties of 
the resulting p-wave supcrfluid below T c . This can be 
done via standard probes of superfluidity, such as the 
presence of vortices in a rotating cloud [29]. Following 
general arguments [30-32] , we expect a p x + ip y ground 
state, i.e. A-f(k) = A Yi,i(k). The anisotropic gapping 
of the spin-j' Fermi surface should yield a signature in 
radio-frequency (RF) spectroscopy, which measures the 
atom transfer rate of one spin species from the interact- 
ing system to an unoccupied energy level [33] , probing the 
spectral function. However, we find that the associated 
peak position in the RF line-shape is at uj ~ Ao(Ao/eFt)j 
a very small energy scale given the smallness of T c com- 
puted above. A more promising route, that we leave for 
future research, is the question of how the onset of pairing 
impacts the formation of the spin- 1 polarons (as reflected 
in, e.g., the spin- 1 RF spectra [8]). 

In summary, we have calculated the induced inter- 
action between like atoms in the normal state of an 
imbalanced two-component Fermi gas. In the absence 
of any competing instabilities (which certainly occur at 
smaller P, where the regimes of magnetic superfluid- 
ity [34], phase separation and, possibly FFLO phase oc- 
cur), this interaction leads to the formation of a p-wave 
supcrfluid in the majority spin species, with a transition 
temperature that peaks, for P close to unity, at a few 
percent of the spin-f Fermi energy. 

This work has been supported by the Louisiana Board 
of Regents, under grant LEQSF (2008-ll)-RD-A-10. 
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